Radio-frequency spectroscopy of weakly bound molecules in spin-orbit coupled atomic 

Fermi gases 
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We investigate theoretically radio-frequency spectroscopy of weakly bound molecules in an ul- 
tracold spin-orbit-coupled atomic Fermi gas. We consider two cases with either equal Rashba and 
Dresselhaus coupling or pure Rashba coupling. The former system has been realized very recently 
at Shanxi University [Wang et ai, arXiv: 1204.1887] and MIT [Cheuk et al, larXiv: 1205.3483] . We 
predict realistic radio-frequency signals for revealing the unique properties of anisotropic molecules 
formed by spin-orbit coupling. 
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I. INTRODUCTION 

The coupling between the spin of electrons to their 
orbital motion, the so-called spin-orbit coupling, lies at 
the heart of a variety of intriguing phenomena in di- 
verse fields of physics. It is responsible for the well- 
known fine structure of atomic spectra in atomic physics, 
as well as the recently discovered topological state of 
matter in solid-state physics, such as topological insu- 
lators and superconductors P, For electrons, spin- 
orbit coupling is a relativistic effect and in general not 
strong. Most recently, in a milestone experiment at the 
National Institute of Standards and Technology (NIST), 
synthetic spin-orbit coupling was created and detected 
in an atomic Bose-Einstein condensate (BEC) of 87 Rb 
atoms [![. Using the same experimental technique, non- 
interacting spin-orbit-coupled Fermi gases of 40 K atoms 
and 6 Li atoms have also been realized, respectively, at 
Shanxi University B and at Massachusetts Institute of 
Technology (MIT) [5!j . These experiments have paved an 
entirely new way to investigate the celebrated effects of 
spin-orbit coupling. 

Owing to the high-controllability of ultracold atoms in 
atomic species, interactions, confining geometry and pu- 
rity, the advantage of using synthetic spin-orbit coupling 
is apparent: (i) The strength of spin-orbit coupling be- 
tween ultracold atoms can be made very strong, much 
stronger than that in solids; (ii) New bosonic topological 
states that have no analogy in solid-state systems may 
be created; (iii) Ultracold atoms are able to realize topo- 
logical superfluids that are yet to be observed in the solid 
state; (iv) Strongly correlated topological states can be 
readily realized, whose understanding remains a grand 
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challenge. At the moment, there has been a flood of the- 
oretical work on synthetic spin-orbit coupling in BECs 
and atomic Fermi gases [l^ - [3ll |. addressing par- 
ticularly new exotic sup erfluid phases arising from spin- 
orbit coupling @, i, Gj| [H . 

In this paper, we investigate theoretically momentum- 
resolved radio-frequency (rf) spectroscopy of an inter- 
acting two-component atomic Fermi gas with spin-orbit 
coupling. The interatomic interactions can be easily ma- 
nipulated using a Feshbach resonance in 40 K or 6 Li atoms 
|32j. It is known that weakly bound molecules with 
anisotropic mass and anisotropic wave function (in mo- 
mentum space) may be formed due to spin-orbit coupling 
[l9l.[23.[23|. Here, we aim to predict observable rf signals 
of these anisotropic molecules. Our calculation is based 
on the Fermi's golden rule for the bound-free transition of 
a stationary molecule [33| . We consider two kinds of spin- 
orbit coupling: (1) the equal Rashba and Dresselhaus 
coupling, Xk x a y , which has been realized experimentally 
at Shanxi University [3] and MIT @, and (2) the pure 
Rashba coupling, \{k y a x — k x a y ), which is yet to be real- 
ized. Here a x and a v are the Pauli matrices, k x and k y 
are momenta and A is the strength of spin-orbit coupling. 
The latter case with pure Rashba spin-orbit coupling is 
of particular theoretical interest, since molecules induced 
by spin-orbit coupling exist even for a negative s-wave 
scattering length above Feshbach resonances [l|l [22], [23[ . 

The paper is structured as follows. In the next section, 
we discuss briefly the rf spectroscopy and the Fermi's 
golden rule for the calculation of rf transfer strength. 
In Sec. Ill, we consider the experimental case of equal 
Rashba and Dresselhaus spin-orbit coupling. We intro- 
duce first the model Hamiltonian and the single-particle 
and two-particle wave functions. We then derive an 
analytic expression for momentum- resolved rf transfer 
strength following the Fermi's golden rule. It can be 
written explicitly in terms of the two-body wave function. 
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We discuss in detail the distinct features of momentum- 
resolved rf spectroscopy in the presence of spin-orbit cou- 
pling, with the use of realistic experimental parameters. 
In Sec. IV, we consider an atomic Fermi gas with pure 
Rashba spin-orbit coupling, a system anticipated to be 
realized in the near future. Finally, in Sec. V, we con- 
clude and make some final remarks. 



of final states for rf transition. Hereafter, without any 
confusion we shall ignore the energy splitting in the bare 
atomic hyperfine levels and set uj^i — 0. To calculate the 
Franck-Condon factor Eq. @, it is crucial to understand 
the initial two-particle bound state | < J?2S ) and the final 
two-particle state |$/). 



II. RADIO-FREQUENCY SPECTROSCOPY 
AND THE FERMI'S GOLDEN RULE 



Radio-frequency spectroscopy, including momentum- 
resolved rf-spectroscopy, is a powerful tool to characterize 
interacting many-body systems. It has been widely used 
to study fermionic pairing in a two-component atomic 
Fermi gas near Feshbach resonances when it crosses from 
a Bardeen-Cooper-Schrieffer (BCS) superfluid of weakly 
interacting Cooper pairs into a BEC of tightly bound 
molecules |34l - l36| . Most recently, it has also been used 
to detect new quasiparticles known as repulsive polarons 
[3^. l38l| . which occur when "impurity" fermionic particles 
interact repulsively with a fermionic environment. 

The underlying mechanics of rf-spectroscopy is sim- 
ple. For an atomic Fermi gas with two hyperfine states, 
denoted as |1) = \f) and |2) = the rf field drives 
transitions between one of the hyperfine states (i.e., \D) 
and an empty hyperfine state |3) which lies above it by 
an energy fno^ due to the magnetic field splitting in bare 
atomic hyperfine levels. The Hamiltonian for rf-coupling 
may be written as, 

V rf = V [ dr y (r) ^ (r) + (r) V>3 (r)l , (1) 



where ipt ( r ) (V4 ( r )) is the field operator which creates an 
atom in |3) (\D) at the position r and, Vb is the strength 
of the rf drive and is related to the Rabi frequency lor 
with Vq = Hujr/2. 

For the rf-spectroscopy of weakly bound molecules that 
is of interest in this work, a molecule is initially at rest in 
the bound state |<I>2_b) with energy Eq = —€b- Here es is 
the binding energy of the molecules. A radio-frequency 
photon with energy ftw will break the molecule and trans- 
fer one of the atoms to the third state |3). In the case 
that there is no interaction between the state |3) and the 
spin-up and spin-down states, the final state f) involves 
a free atom in the third state and a remaining atom in 
the system. According to the Fermi's golden rule, the rf 
strength of breaking molecules and transferring atoms is 
proportional to the Franck-Condon factor [33j . 



F{u) = \{* f \ V rf \$ 2B )\ 2 6 



U3± 



E f - E 



(2) 



where the Dirac delta function ensures energy conserva- 
tion and Ef is the energy of the final state. The inte- 
grated Franck-Condon factor over frequency should be 
unity, f^°°F(uj)dijj = 1, if we can find a complete set 



III. EQUAL RASHBA AND DRESSELHAUS 
SPIN-ORBIT COUPLING 



Let us first consider a spin-orbit-coupled atomic Fermi 
gas realized recently at Shanxi University |4| and at MIT 
[5j. In these two experiments, the spin-orbit coupling is 
induced by the spatial dependence of two counter propa- 
gating Raman laser beams that couple the two spin states 
of the system. Near Feshbach resonances, the system may 
be described by a model Hamiltonian % = %o + Hint, 
where 



Ho = E/ dr ^W 



h 2 k 2 

2m 



•Vv (r) 



dr 



,frt (p) (fke«**») 



i> i (r) + H.c. 



(3) 



is the single-particle Hamiltonian and 



H int = U dr^\ (r) ( r ) ^ ( r ) ^ ( r) ( 4 ) 



is the interaction Hamiltonian describing the contact in- 
teraction between two spin states. Here, (r) is the 
creation field operator for atoms in the spin-state a, 
hk = —iKV is the momentum operator, fiij is the cou- 
pling strength of Raman beams, • = 2tt /\r is deter- 
mined by the wave length \r of two Raman lasers and 
therefore 2hkn is the momentum transfer during the two- 
photon Raman process. The interaction strength is de- 
noted by the bare interaction parameter U$. It should be 
regularized in terms of the s-wave scattering length a s , 
i.e., 1/U = m/(4nh 2 a s ) - £ k m/(h 2 k 2 ). 

To solve the many-body Hamiltonian Eq. ©, it is 
useful to remove the spatial dependence of the Raman 
coupling term, by introducing the following new field op- 
erators %jj a via 



^ t (r) = e +ikRX ^ (r) 

-iknx n 



^(r) = e-^^(r). 



(5) 
(G) 



With the new field operators ip a , the single-particle 
Hamiltonian then becomes. 



(k ± fc_Re 3 
2m 



-4>o- (r) 



4>\ (r) $± (r) + H.c. 



(7) 
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where in the first term on the right hand side of the 
equation we take "+" for spin- up atoms and "— " for spin- 
down atoms. The form of the interaction Hamiltonian is 
invariant, 



where ^| (r) = E q c q3 eiqr and fa ( r ) = E q c qi eiq r - 
Hereafter, we shall denote Ck3 and Cko- as the field op- 
erators (in the momentum space) for atoms in the third 
state and in the spin-state a, respectively. 



H in t = U / dr $ (r) (r) fa (r) fa (r) 



(8) 



However, the rf Hamiltonian acquires an effective mo- 
mentum transfer kne x , 



V rf = V / dr \e~ ikRX fa 3 (r) fa (r) + H.c. 



(9) 



For later reference, we shall rewrite the rf Hamiltonian 
in terms of field operators in the momentum space, 



Vrf = VoJ2 ( C q-fcfle x ,3 C qi + H ' C 



(10) 



A. Single-particle solution 



Using the Pauli matrices, the single-particle Hamilto- 
nian takes the form. 



Ho= / di#I(r),^(r)] 



h 2 (fc| + k 2 
2m 



+ ha x + Xk x a 2 



V't ( r ) 

L fa to 



(ii) 



where for convenience we have defined the spin-orbit coupling constant denoted as A = h 2 kn/m and an "effective" 
Zeeman field h = Qr/2. This Hamiltonian is equivalent to the one with equal Rashba and Dresselhaus spin-orbit 
coupling, Xk x a y . To see this, let us take the second transformation and introduce new field operators (r) via 



^t(r) = _[$ t ( r )-<^(r)], 

fair) = ^|[*t(r)+^(r)]. 
Using these new field operators, the single-particle Hamiltonian now takes the form, 



Ho = J dr[*\ (r) , *j (r) ] 



fi 2 (fc| + fc 2 ) 



2m 



\k x (j y + ha z 



*t ( r ) 
*l(r) 



(12) 
(13) 

(14) 



which is precisely the Hamiltonian with equal Rashba 
and Dresselhaus spin-orbit coupling. 

The single-particle Hamiltonian Eq. (fTTj) can be diag- 
onalized to yield two eigenvalues 

Here "±" stands for the two helicity branches. The single- 
particle eigenstates or the field operators in the helicity 
basis take the form 

Ck+ = +c kt cos6> k + CkiSin6» k , (16) 

Ck- = -Cktsin^k + Ck^cos^k, (17) 

where 

<9 k = arctan[(v//i 2 + A 2 fc 2 - Xk x )/h] > (18) 



is an angle determined by h and k x . Note that, 

cos 2 k = if 1 + ; Xkx 1 , (19) 
sin 2 ^ = 1(1- - Xkx 1 . (20) 

Note also that the minimum energy of the single-particle 
energy dispersion is given by [24j 

_ h 2 k\ m\ 2 h 2 h 2 h 2 h 2 

min ~ 2m 2h 2 2mA 2 ~ 2mA 2 ' ( ' 
if h < mX 2 /h 2 . 
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B. The initial two-particle bound state |$2b) singlet and triplet components fl9i . |22| . |23] . The wave 

function at zero center-of-mass momentum, |$2s)j may 
In the presence of spin-orbit coupling, the wave func- be written as [22j|, 
tion of the initial two-body bound state has both spin 



*2s) = [^U ( k ) 4t c -l4 + ( k ) c L c -kt + V^tt ( k ) 4t c -kt + ( k ) c L c -kJ l vac ) 



(22) 



where and cj^ are creation field operators of spin-up and spin-down atoms with momentum k and C = 

Skll^Ti ( k )| 2 + l' ! /'lt( k )| 2 + IV'ttl' 4 )! 2 + ( k )| 2 ] i s the normalization factor. From the Schrodinger equation 
{Ho+Hint) | < J ) 2_b) = Eo \$2b), we can straightforwardly derive the following equations for coefficients tp aa > appearing 
in the above two-body wave function [22] : 



E 



h 2 k 2 R 



h 2 k 2 R 



Eq 

Eq 



h 2 k 2 

m 

h 2 k 2 

m 
h 2 k 2 R 



h 2 k% 

m 



2Xk x 

2Xk x 

H 2 k 2 
m 

h 2 k 2 

m 



i> n (k) = + y E &n ( k ') - ( k ')l + ^tt ( k ) + W u (k) , 

k' 

^t(k) = -^^[^ n (k')-^t(k')] + ^tt( k ) + ^u( k )» 

k' 

t/j n (k) = hip n (k) + hip lf (k) , 

i> u (k) = hip n (k) + hijftf (k) , 

I 



(23) 

(24) 

(25) 
(26) 



where Eq = —£b < is the energy of the two-body 
bound state. Let us introduce A k = —£b — (h 2 k R /m + 
h 2 k 2 /m) < and different spin components of the wave- 
functions, 



1>* ( k ) 
1>a ( k ) 
It is easy to see that, 

Xk 



-±= (k) - Vit ( k )l , (27) 
-i=[^ ti (k) + Vit( k )]- (28) 



and 



ip n (k) 

V>U ( k ) 
% ( k ) 

-4, 



1>a (k) , 

i>a (k) , 



Au - 2h Ai, 



2h 



(29) 
(30) 

V^( k ), (31) 



AX 2 k 2 



A k - Ah 2 /A k 



^(k) = C/ ]T^( k ')- (32) 



As required by the symmetry of fermionic system, the 
spin-singlet wave function tp s (k) is an even function 
of the momentum k, i.e., ip s (— k) = ?/> s (k), and the 
spin-triplet wave functions are odd functions, satisfy- 
ing Va(-k) = -1p a (k), ^tt( _k ) = _ ^tt( k ) aild 



4>ll (~ k ) = — ^44 ( k )- The un-normalized wavefunction 
</> s (k) = [A k - 4A 2 fc 2 /(A k - 4^/Ak)]- 1 is given by, 



1>s ( k ) 



1 



/i 2 + A 2 fc 2 



A k 



X 2 k 2 x A w 



A 2 



(33) 



4(h 2 + X 2 k 2 )_ 

Using Eq. (|32|) and un-normalized wave function ip s (k), 
the bound-state energy Eq or the binding energy es is 
determined by Uq ^ k ^ s (k) = 1, or more explicitly, 



rn 



Airh 2 



0. 



(34) 



Here we have replaced the bare interaction strength Uq 
by the s-wave scattering length a s using the standard 
regularization scheme mentioned earlier. The normaliza- 
tion factor of the total two-body wave function is given 
by, 



c = x> a (k)r 



i 



2A 2 fc 2 



2A 2 fc 2 



{A k - 2hf {A k + 2hf 



(35) 



C. The final two-particle state $/) 



Let us consider now the final state |$/). For this pur- 
pose, it is useful to calculate 



V rf |#2B> = V E ^q-fe^^-qi I*™) 



(36) 
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and then determine possible final states. It can be readily 
seen that, 



Vr f |$ 2 B> = - - |L E c iq-^, 3 {&u (q) - (-q)] 4 + (q) - ^ (-q)] 4} 



(37) 



Rewriting and ipy* in terms of -0 S and tp a as shown in Eqs. (|27[) and (|28[) . and exploiting the parity of the wave 
functions, we obtain a general result valid for any type of spin-orbit coupling, 



v rf i# afl ) = -y^v ^ c+ q _ fej?ex3 {[^ s (q) + ^ a (q)] 4 + v^Vw (q) 4} l vac > 



(38) 



To proceed, we need to rewrite the field operators 4- and 4, in terms of creation operators in the helicity basis. For 
the case of equal Rashba and Dresselhaus spin-orbit coupling, using Eqs. (TT5|) and (fTT|) , we find that 



Thus, we obtain 



cos V q c qH 
y q c q4 



Vrf\^B)^-J ? V J2^ kR 



sin6» q c q _, 
cos 6» q c q _ . 



(39) 
(40) 



(41) 



where 



Sq+ = h/'s (q) + -0a (q)] cos 6q + V2ipa (q) sin6> q 
s q - = (q) + V'a (q)] sin 6> q - (q) cos 6> q 



(42) 
(43) 



Eq. (|4ip can be interpreted as follows. The rf photon breaks a stationary molecule and transfers a spin-down atom to 
the third state. We have two possible final states: (1) we may have two atoms in the third state and the upper helicity 
state, respectively, with a possibility of |s q+ | /C; and (2) we may also have a possibility of |s q _| jC for having two 
atoms in the third state and the lower helicity state, respectively. 



D. Momentum-resolved rf spectroscopy 



Taking into account these two final states and using the Fermi's golden rule, we end up with the following expression 
for the Franck-Condon factor, 



s q+ <5 1 w 



i , 2 r ( £q- 

1 + S q _(3 | LO 



(44) 



where 



(4 



2 m 



± V^ 2 + A 2 ^ 



fr 2 (q + fegga)^ 
2m 



(45) 



The two Dirac delta functions in Eq. (|44[) are due to energy conservation. For example, the energy of the initial state 
(of the stationary molecule) is Eq = —eb, while the energy of the final state is ?i 2 (q + kjie x ) 2 / (2m) for the free atom 



in the third state and h (fcjj. 



7 2 )/(2m) 



A 2 <7 2 for the remaining atom in the upper branch. Therefore, the 



rf energy Huj required to have such a transfer is given by £ q +, as shown by the first Dirac delta function. It is easy to 

check that the Franck-Condon factor is integrated to unity, J_°° F (ui) = 1. 

Experimentally, in addition to measuring the total number of atoms transferred to the third state, which is propor- 
tional to F(uj), we may also resolve the transferred number of atoms for a given momentum or wave-vector k x . Such 
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a momentum-resolved rf spectroscopy has already been implemented for a non-interacting spin-orbit coupled atomic 
Fermi gas at Shanxi University and at MIT. Accordingly, we may define a momentum-resolved Franck-Condon factor, 



F(k x ,Lj) = 



-y 

C 



q+ 



Sq_S I LJ - 



e, 



(46) 



where the summation are now over the wave- vector = (q y , q z ) and we have defined k x = q x + Icr by shifting the 
wave- vector q x by an amount k R . This shift is due to the gauge transformation used in Eqs. ([5]) and ([5]). With the 
help of the two Dirac delta functions, the summation over qj_ may be done analytically. We finally arrive at, 



F(k x ,oj) 

where (x) is the step function and 



8n 2 hC 



K + e ( q l + ) + sl_e 



q\. ± = j(lJ - H) - (kl + ql + q x k R ± ™JWTWl 



(47) 



(48) 



It is understood that we will use q = (q x ,q± t +) in the 
calculation of s q+ and q = (q x , q± —) in s q _. 

We may immediately realize from the above expres- 
sion that the momentum-resolved Franck-Condon factor 
is an asymmetric function of k x , due to the coexistence of 
spin-singlet and spin-triplet wave functions in the initial 
two-body bound state. Moreover, the contribution from 
two final states or two branches should manifest them- 
selves in the different frequency domain in rf spectra. As 
we shall see below, these features give us clear signals of 
anisotropic bound molecules formed by spin-orbit cou- 
pling. On the other hand, from Eq. (|4T|) . it is readily 
seen that once the momentum-resolved rf spectroscopy 
is measured with high resolution, it is possible to deter- 



mine precisely s'i, and 



and then re-construct the 



two-body wave function of spin-orbit bound molecules. 



E. Numerical results and discussions 

For equal Rashba and Dresselhaus spin-orbit coupling, 
the bound molecular state exists only on the BEC side 
of Feshbach resonances with a positive s-wave scattering 
length, a s > [24]. Thus, it is convenient to take the 
characteristic binding energy Eb = fi 2 / (ma 2 ) as the unit 
for energy and frequency. For wave- vector, we use kn — 
m\/h 2 as the unit. The strength of spin-orbit coupling 
may be measured by the ratio 



Ex 
E B 



mAa, 



(49) 



where we have defined the characteristic spin-orbit en- 
ergy E\ = m\ 2 /h 2 = f?k\Jm. Note that, the spin-orbit 
coupling is also controlled by the effective Zeeman field 
h = Qfj/2. In particular, in the limit of zero Zeeman field 
il/f = 0, there is no spin-orbit coupling term as shown 
in the original Hamiltonian Eq. ([3|), although the char- 
acteristic spin-orbit energy E\ ^ 0. Using kf, and Eb 
as the units for wave-vector and energy, we can write 



a set of dimensionless equations for the binding energy 
eb = —Eq, normalization factor C, Franck-Condon factor 
F (oj) and the momentum-resolved Franck-Condon fac- 
tor F (u!,k x ). We then solve them for given parameters 
E\/Eb and h/ E\. In accord with the normalization con- 
dition F (ui) = 1, the units for F (lj) and F (k x ,ui) 
are taken to be 1/Eb and I/^s/cr), respectively. 



0.3 



3 

ST 



,0.2 



0.1 



0.0 



EJE =2 




<o/E„ 



Figure 1: (color online) Franck-Condon factor of weakly 
bound molecules formed by equal Rashba and Dresselhaus 
spin-orbit coupling, in units of Eg 1 . Here we take h — Ex/2 
or Qr = h 2 k R /m and set E X /E B = 0.5, 1, and 2. The re- 
sult without spin-orbit coupling is plotted by the thin dashed 
line. The Inset shows the different contribution from the two 
final states at E\/Eb = 1. The one with a remaining atom 
in the lower (upper) helicity branch is plotted by the dashed 
(dot-dashed) line. 

Fig. I displays the Franck-Condon factor as a func- 
tion of the rf frequency at h/E\ — 0.5 and at several 
ratios of E\/Eb as indicated. For comparison, we show 
also the rf line-shape without spin-orbit coupling (33| . 
F (lj) — (2/tt)\/lj — Eb/lj 2 , by the thin dashed line. In 
the presence of spin-orbit coupling, the existence of two 
possible final states is clearly revealed by the two peaks 
in the rf spectra. This is highlighted in the inset for 
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E\/Eb — 1, where the contribution from the two possi- 
ble final states is plotted separately. The main rf response 
is from the final state with the remaining atom staying 
in the lower helicity branch, i.e., the second term in the 
Franck-Condon factor Eq. (pPfj) . The two peak positions 
may be roughly estimated from Eq. (j48|) for the threshold 
frequency lo c ± of two branches, 



tko, 



c± 



ft 2 (k% + q 2 + q x k R ) 



rn 



±^h 2 Tx 2 q 2 



With increasing spin-orbit coupling, the low- frequency 
peak becomes more and more pronounced and shifts 
slightly towards lower energy. In contrast, the high- 
frequency peak has a rapid blue-shift. 




" 3 1 

1 4 
o4 



Figure 2: (color online) Linear contour plot of momentum- 
resolved Franck-Condon factor, in units of (-Bsfcii)" 1 . Here 
we take h = E\/2 and consider E\/Eb = 0.1, 0.5, 1, and 2. 




Figure 3: (color online) Energy distribution curve of 
the momentum-resolved Franck-Condon factor, in units of 
(-Esfcij) -1 . We consider several values of the rf frequency w as 
indicated, under given parameters h = E\/2 and E\ = Eb- 

Fig. 2 presents the corresponding momentum-resolved 
Franck-Condon factor. We find a strong asymmetric dis- 



tribution as a function of the momentum k x . In par- 
ticular, the contribution from two final states are well 
separated in different frequency domains and therefore 
should be easily observed experimentally. The asymmet- 
ric distribution of F(k Xl cu) is mostly evident in energy 
distribution curve, as shown in Fig. 3, where we plot 
F (k x ,uj) as a function of k x at several given frequencies 
uj. In the experiment, each of these energy distribution 
curves can be obtained by a single-shot measurement. 




Figure 4: (color online) Zeeman-field dependence of the 
Franck-Condon factor at E\ = Eb- Here we vary the effective 
Zeeman fields, h/E\ = 0.4, 0.6, 0.8 and 1.0. The inset shows 
the momentum-resolved Franck-Condon factor at h = E\. 

We finally discuss the effect of the effective Zeeman 
field h = CIr/2. Fig. 4 shows how the line-shape of 
Franck-Condon factor evolves as a function of the Zee- 
man field at E\/Eb = 1. In general, the larger Zeeman 
field the stronger spin-orbit coupling. Therefore, as the 
same as in Fig. 1, the increase in Zeeman field leads to 
a pronounced peak at about the binding energy. There 
is a red-shift in the peak position as the binding energy 
becomes smaller as the Zeeman field increases. As antic- 
ipated, the larger the Zeeman field, the more asymmetric 
F{k x ,uj) becomes. In the inset, we show as an example 
the contour plot of F (k x ,u>) at h/E\ = 1. 



IV. RASHBA SPIN-ORBIT COUPLING 

We now turn to the case with pure Rashba spin-orbit 
coupling, X(k y (T x — k x a y ), which may be realized exper- 
imentally in the near future. The single-particle Hamil- 
tonian may be written as [23|, 



Hq 



dr 



4(r),Vj (r) 



S 



V> t (r) 
^ (r) 



where the matrix 
S 



k 2 ) I (2m) 



h 2 (k 2 R 



l\ iky ^) 

k 2 ) I (2m) 



(51) 



(52) 



Here A is the coupling strength of Rashba spin-orbit cou- 
pling, kji = mX/h 2 , and we have added a constant term 



s 



h k R /(2m) to make the minimum single-particle energy Here ip^ = SLig(k x , k y ) is the azimuthal angle of the wave- 



zero |24j, i.e., E min = 0. 



vector kj_ in the x — y plane. 



A. Single-particle solution 

We diagonalize the matrix S to obtain two helicity 
eigenvalues [23^ . 



h 2 (kj, + k 2 ) 

Am 



(53) 



where k± = yk 2 + k 2 and "±" stands for the two helicity 

branches. For later reference, the field operators in the 
original spin basis and in the helicity basis are related by, 



'kf 



V2 



<■;., = — = 1 4 , + • 



i<Pk r f 

k- I ' 



'k, = ^ (--"^i , + <{- 



(54) 
(55) 



B. The initial two-particle bound state |4?2b) 



In the case of Rashba spin-orbit coupling, the two-body 
wave function can still be written in the same form as 
in the previous case, i.e., Eq. (f2"2")l . and the Schrodinger 
equation leads to [22j, 



A^ u (k) = +^ E &n (k') - (k')] - A {ky - ik x ) ^ t (k) + A (ky + ik x ) j, u (k) , 



U, 



A^t ( k ) - E &n (k') - ^ (k')] + A (k y - ik x ) </> n (k) - A (ky + ik x ) V u (k) , 



A k ip n (k) = -\(k y + ik x )ip n (k) + X(k y 
Ak^u ( k ) = + A ( k y - ik x) i>n (k) - A (k y 



■ik x )t/)tf (k), 
ik x )iptf ( k ), 



(56) 

(57) 

(58) 
(59) 



r 



where A k = Eq — (t?k 2 R jm + h 2 k 2 /m) < 0. It is easy to 
show that ip a (k) = and 



-4, 



Thus, we obtain the (un-normalized) wavefunction: 



i>* (k) 



1 



1 



and the equation for the energy Eq, 

1/2 1/2 



(61) 



E 



So - 2E U+ Eq - 2£ k - h 2 k 2 



The spin-triplet wave functions if?^ (k) and ipn (k) are 
given by, 



ip n (k) 
V-U (k) 



V2Xk± 
Eq — 2ek 

V2\k± 

Eq — 2fik 



^ (k) , (63) 
V> s (k) , (64) 



where £k = fo 2 k 2 /(2m). The normalization factor for the 
two-body wave function is therefore, 



c = £>. (k)l 2 



1 + 



A\ 2 k 



2 1.2 



(E - 2e k r 



(65) 



C. The flnal two-particle state $/) 

To obtain the final state, we consider again V r f | < &2s) ■ 
In the present case, we assume that the rf Hamiltonian 
is given by, 



\>rf = Vq E ( C q3 C qi + C qi C 1 3 



(66) 



Following the same procedure as in the case of equal 
Rashba and Dresselhaus coupling, it is straightforward 
to show that, 
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V rf |*2B) = -\j^Y. °iq3 h (q) C qt + (q) <ij I vac) 

q 

Using Eqs. f)54|) and (|55p to rewrite cLj. and in terms of and Cq_, we obtain, 



(67) 



q L 



-q3 L q 



q3 c q- 



Eq — 2£?q_|_ _Eo — 2Sq_ 



vac) 



(68) 



Therefore, we have again two final states, differing in the helicity branch that the remaining atom stays. The 
remaining atom stays in the upper branch with probability (2C)~ 1 (Eq — 2_E q+ )~ 2 , and in the lower branch with 
probability (2C)- 1 (E - 2E^)- 2 . 



D. Momentum-resolved rf spectroscopy 

Using the Fermi's golden rule, we have immediately the Franck-Condon factor, 



2(e B +2E k+ f 2(e B + 2E^) 2 



(69) 



where 



h 2 k% h 2 k 2 ^ . 

£ k± = e B + + ± Xk ± . 



2m m 

For Rashba spin-orbit coupling, it is reasonable to define the following momentum-resolved Franck-Condon factor 

S (uj — £k+ /ft) 6 (oj — £k- / K) 



(70) 



F(k ± ,u) = ±J2 



2(e B + 2E k+ ) 2 2{e B + 2E^f 



(71) 



where we have summed over the momentum k z . Integrating over k z with the help of the two Dirac delta functions, 
we find that, 



F(k±,u) 



16n 3 hC 



Q(k 2 z+ ) 



Q(kl_) 



(Hlu + h 2 k 2 R /2m + \k±) k z+ (hid + h 2 k 2 R /2m - Xk±) k z 



where 



2 m ( e B 

z± = ~ft\ w ~~ft 



^f + kl± k R k ± 



(72) 



(73) 



It is easy to see that the threshold frequencies for the two 
final states are given by, 



huj c+ 

hUJr- 



(-B + 



ft 2 k\ 
2m ' 

H 2 k 2 R 
Am ' 



(74) 
(75) 



which differ by an amount of ft 2 k R /(4m) = Ex/A. Near 
w c _, we find approximately that F (w) oc 0(w — w c _)/w 2 . 
Thus, the lineshape near the threshold is similar to that 
of a two-dimensional (2D) Ferm gas [36J. This similarity 
is related to the fact that at low energy a 3D Fermi gas 



with Rashba spin-orbit coupling has exactly the same 
density of states as a 2D Fermi gas [lij . 



E. Numerical results and discussions 

For the pure Rashba spin-orbit coupling, the molec- 
ular bound state exists for arbitrary s-wave scattering 

length a s [H, [HHH]. We sha11 take k R = ™V^ 2 and 
E\ = m\ 2 /h 2 as the units for wave- vector and energy, 
respectively. With these units, the dimensionless inter- 
action strength is given by h 2 /(mXa s ). The spin-orbit 
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effect should be mostly significant on the BCS side with 
H 2 1 (mXa s ) < 0, where the bound state cannot exist with- 
out spin-orbit coupling. 




Figure 5: (color online) Franck-Condon factor of weakly 
bound molecules formed by Rashba spin-orbit coupling, in 
units of -E^" 1 . Here we take h 2 / \m\a s ) = — 1 (dashed line), 
(solid line), and 1 (dot-dashed line). In the deep BCS limit, 
ti 2 1 (mXa B ) — > — oo, the Franck-Condon factor peaks sharply 
at frw ~ E\/2 and becomes a delta-like distribution. 

Fig. 5 shows the Franck-Condon factor at three dif- 
ferent interaction strengths H 2 / (mXa s ) = —1, 0, and +1. 
The strong response in the BCS regime (a s < 0) or in 
the unitary limit (a s — > ±oo) is an unambiguous signal 
of the existence of Rashba molecules. In particular, the 
rf line-shape in the BCS regime shows a sharp peak at 
about Huj ~ E\/2 and decays very fast at high frequency. 

In Fig. 6, we present the corresponding momentum- 
resolved Franck-Condon factor F (k±,u), in the form of 
contour plots. We can see clearly the different response 
from the two final states. The momentum-resolved rf 
spectroscopy is particularly useful to identify the contri- 
bution from the final state that has a remaining atom in 
the upper branch, which, being integrated over fcj_, be- 
comes too weak to be resolved in the total rf spectroscopy. 
Finally, we report in Fig. 7 energy distribution curves of 
F(k±,u) in the unitary limit H 2 /(mXa s ) = 0. We find 
two sharp peaks in each energy distribution curve, aris- 
ing from the two final states. When measured experi- 
mentally, these sharp peaks would become much broader 
owing to the finite experimental energy resolution. 



V. CONCLUSIONS 

In conclusions, we have investigated theoretically the 
radio-frequency spectroscopy of weakly bound molecules 
in a spin-orbit coupled atomic Fermi gas. The wave 
function of these molecules is greatly affected by spin- 
orbit coupling and has both spin-singlet and spin-triplet 
components. As a result, the line-shape of the to- 
tal radio-frequency spectroscopy is qualitatively differ- 
ent from that of the conventional molecules at the BEC- 




0.001 



klk 



Figure 6: (color online) Contour plot of momentum-resolved 
Franck-Condon factor of weakly bound molecules formed by 
Rashba spin-orbit coupling, in units of (Exkn)~ . The in- 
tensity increases from blue to red in a logarithmic scale. We 
consider h 2 / (m\a s ) = — 1 (a), (b), and 1 (c). 



BCS crossover without spin-orbit coupling. In addi- 
tion, the momentum-resolved radio-frequency becomes 
highly asymmetric as a function of the momentum. 
These features are easily observable in current experi- 
ments with spin-orbit coupled Fermi gases of 40 K atoms 
and 6 Li atoms. On the other hand, from the high- 
resolution momentum-resolved radio-frequency, we may 
re-construct the two-body wave function of the bound 
molecules. 

We consider so far the molecular response in the radio- 
frequency spectroscopy. Our results should be quantia- 
tively reliable in the deep BEC limit with negligible num- 
ber of atoms, i.e., in the interaction parameter regime 
with l/(kpa s ) > 2. However, in a real experiment, in or- 
der to maximize the spin-orbit effect, it is better to work 
closer to Feshbach resonances, i.e., l/(kpa s ) ~ 0.5. Un- 
der this situation, the spin-orbit coupled Fermi gas con- 
sists of both atoms and weakly bound molecules, which 
may strongly interact with each other. Our prediction 
for the molecular response is still qualitatively valid, with 
the understanding that there would be an additional pro- 
nounced atomic response in the rf spectra. A more in- 
depth investigation of radio-frequency spectroscopy re- 
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Figure 7: (color online) Momentum-resolved Franck-Condon 
factor of Rashba molecules at h 2 / (m\a s ) = 0, shown in the 
form of energy distribution curves at several rf frequencies as 
indicated. 



quires complicated many-body calculations beyond our 
simple two-body picture pursued in the present work. 
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